Local and chain dynamics in miscible polymer blends: A Monte Carlo simulation 

study 
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Local chain structure and local environment play an important role in the dynamics of polymer 
chains in miscible blends. In general, the friction coefficients that describe the segmental dynamics 
of the two components in a blend differ from each other and from those of the pure melts. In 
this work, we investigate polymer blend dynamics with Monte Carlo simulations of a generalized 
bond-fluctuation model, where differences in the interaction energies between non-bonded nearest 
neighbors distinguish the two components of a blend. Simulations employing only local moves and 
respecting a non-bond crossing condition were carried out for blends with a range of compositions, 
densities, and chain lengths. The blends investigated here have long-chain dynamics in the crossover 
region between Rouse and entangled behavior. In order to investigate the scaling of the self-diffusion 
coefficients, characteristic chain lengths iV c are calculated from the packing length of the chains. 
These are combined with a local mobility fj, determined from the acceptance rate and the effective 
bond length to yield characteristic self-diffusion coefficients D c = fi/N c . We find that the data for 
both melts and blends collapse onto a common line in a graph of reduced diffusion coefficients D/D c 
as a function of reduced chain length N/N c . The composition dependence of dynamic properties is 
investigated in detail for melts and blends with chains of length twenty at three different densities. 
For these blends, we calculate friction coefficients from the local mobilities and consider their com- 
position and pressure dependence. The friction coefficients determined in this way show many of 
the characteristics observed in experiments on miscible blends. 
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I. INTRODUCTION 

Processes on different length scales affect dynamic 
properties of polymer melts and blendsi2i2i4 The dy- 
namics of miscible blends have been the subject of a large 
number of recent in vesti gations with experimental (see, 
for example, Refs. 

G3IH, simulation (cf. Refs. ll9l2M2]l2a2a232"5l l . and 

theoretical methods (cf. Refs. OsaE 

Monomeric friction coefficients, which are in- 
versely proportional to the mobility of short chain seg- 
ments, are a convenient way to describe the effect of local 
dynamic properties on global dynamic properties such as 
the viscosity^ 

Since blending changes the local environment of 
the chain segments of a polymer it has a strong ef- 
fect on the local dynamics of the chains. From 
experimental 5 ! 6 ! 7 ! 8 ! 9 ! 10 . 11 ! 12 . 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 and simula- 
tion worh 19 ! 20 ! 21 ! 22 ! 23 ! 24 ! 25 it is found that the local dy- 
namics of the two blend components differ from each 
other and the pure melts. The addition of slow (high 
friction coefficient) component to a blend is found to in- 
crease the friction coefficients of both components and, 
conversely, the addition of fast component is found to 
speed up both components. The effects of blending are 
most pronounced near the glass transition, however, they 
are observed even at high temperatures 1 ^ 4 * 1 ^ 9 " and in 
blends where one of the components is dilute i 1 ^ 6 *^ 

Several recently developed models for miscible poly- 



mer blends relate differences in the component dy- 
namics to local variations in the glass transition 
temperature induced by local variations in blend 
composition; 8 i 26 i 27 i 28 i 29 i 30 i 31 i 32 while others consider "in- 
trinsic" effects, due to differences in the chain structure 
of the two components, in addition to local density and 
composition variations^i 3 ^* 3 ^ Unfortunately, it is still 
difficult, in general, to predict the local friction coeffi- 
cients in a blend from those of the melts. 

Self diffusion coefficients give information about the 
chain dynamics of polymer melts and blends. The Rouse 
model prediction for the self-diffusion coefficient of a 
polymer chain may be written asi^ 



(N 



where T is the temperature, fee is Boltzmann's constant, 
N is the chain length, £ is the so-called monomeric fric- 
tion coefficient, and fi = k^T/C, is the corresponding mo- 
bility. The long-time dynamics of long polymer chains are 
dominated by entanglement effects, which are not part of 
the Rouse model. The reptation modeiiiSi 3 . describes the 
entanglements of a chain with other chains in terms of 
a tube, which restricts the motion of the chain perpen- 
dicular to the tube. If the average number of monomers 
between entanglements is denoted by N c , then the repta- 
tion prediction for the self-diffusion coefficient takes the 
form 

k B TN c 1N C 
d g = = -^fD R . (2) 
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The entanglement length of a polymer melt may be de- 
termined directly from experimental data on the plateau 
modulusii^ For viscoelastic properties, it differs by a 
constant prefactor from the characteristic chain length 
N c , that separates short chain (unentangled) behav- 
ior from long chain (entangled) behavior^ Experi- 
mental, theoretical and simulation work on polymer 
melts suggests that the tube diameter of the repta- 
tion theories is proportional to the so-called packing 
length, 36 i 3 Ti 38 i 3fl i 4n i 41 which may be defined as 3 ^ 
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(3) 



where p is the monomer density and R^, is the average 
squared end-to-end vector of the chains. Hence, the en- 
tanglement length is expected to be proportional to the 
chain length corresponding to p, which yields* 



Nr. 
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The transition between the unentangled and entangled 
regimes is not sharp and simulation data are often found 
to be in the crossover region between unentangled and 
reptation behaviori 42 i 43 i 44 i 45 i 46 i 47 In the following, it is 
convenient to rescale the chain length iV by the char- 
acteristic chain length N c cx N e and the self-diffusion 
coefficient I? by a characteristic diffusion coefficient D c 
defined as— 



D c = ±. 

N 



(5) 



In this way, the predictions for the self-diffusion coeffi- 
cient may be summarized as 



D 
D~c 



(N/N c )~ for 7V«iV c 
(N/N c )~ 2 for N~»N C 1 ' 



so that, at the chain length N c , the extrapolations from 
both power laws yield the same value, D/D c — 1. For 
long but not completely entangled chains, Hess— ^ has 
proposed a crossover equation for the self-diffusion coef- 
ficient, which may be written in our notation as 



D_ = (N/N c )~ 
D c 1 + N/N c 



(7) 



This expression has been found to give a good represen- 
tation of data for self-diffusion coefficients of the bond- 
fluctuation model. 50 However, the resulting values for the 
characteristic chain length N c are considerably smaller 
than those obtained by superimposing simulation data 
with experimental data that extends deeply into the en- 
tangled regime— 

Simulation work on polymer blends has been car- 
ried out with atomistic————— 5 - and coarse grained 
models— i— Molecular dynamics simulations of atomistic 
models give access to the chain structure and dynamics 
of realistic polymers and allow a detailed investigation of 



the environments of the chain segments^ 5 Simulations of 
coarse grained models, on the other hand, can be used 
to isolate particular effects such as isotope effects— and 
differences in chain stiffness— 

In this work, we investigate chain and local dynamics 
of miscible polymer blends with the aid of Monte Carlo 
simulations of a lattice model. We are interested in en- 
ergetic effects and represent the two components of a 
blend by chains that differ only in interaction energies. 
The model is a modification of Shaffer's bond-fluctuation 
model for athermal melts 5 ——— and is introduced in Sec- 
tion [H] Section [H] also provides some details about the 
Monte Carlo simulations. In Section Infl we discuss the 
evaluation of the Monte Carlo simulations and describe 
the construction of scaling variables for the self-diffusion 
coefficients. Results of our simulations are presented in 
Section Hvl and discussed in Section [V] 



II. BLEND MODEL FOR MONTE CARLO 
SIMULATION 

Shaffer's bond fluctuation model— ——^ is a lattice 
model for polymer chains, where the monomers occupy 
sites of a simple cubic lattice. In the following, we will 
take the size of the unit cell, i.e. the lattice constant a, 
as the unit for length. The monomers are connected by 
bonds of three possible lengths, namely 1, y/2, and \/3, 
corresponding to the sides, face diagonals, and body di- 
agonals of the unit cell of the lattice. Monte Carlo sim- 
ulations of the model employ only local moves, where an 
attempt is made to displace a randomly chosen monomer 
by one lattice site along any of the three coordinate di- 
rections. One attempted elementary move per monomer 
in the system is called one Monte Carlo step (MCs) and 
will be our unit of time. Shaffer considered two versions 
of this model. In the first version, bonds are allowed to 
cross each other with the result that the chains do not 
entangle; in the second, bond crossings are prohibited 
and entanglement effects become apparent. We adopt 
the second approach and prohibit bond crossings in this 
work. 

Shaffer's model describes athermal melts; the 
monomers interact only through hard core repulsion, 
which is enforced by prohibiting double occupation of 
lattice sites— In this work, we are interested in blends 
of polymers that have identical bond structures but dif- 
fer in their monomer-monomer interactions. To this 
end, we add attractive interactions between non-bonded 
monomers occupying nearest-neighbor sites to the model, 
see Fig.Q] For a binary blend of chains of type A and B, 
the energy parameters €aa , £bb, and €ab , describe inter- 
actions between monomers of type A, type B, and mixed 
interactions, respectively. The total internal energy of 
the system is given by 



E = Naa^aa + Nbb^bb + Nab^ab, 



(8) 



where A£y, i,j € {A, £?}, denotes the number of near- 
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FIG. 1: Illustration of the generalization of Shaffer's bond 
fluctuation model to polymer blends. The circles connected 
by heavy lines represent polymer chains in a blend, gray 
for component A and black for component B. Allowed bond 
lengths are 1, \/2 and \/3 corresponding to the sides, face 
diagonals, and body diagonals of the unit cell of the simple 
cubic lattice. The arrows show examples of nearest neigh- 
bor interactions between like (caa and ess) and unlike (cab) 
monomers. 

est neighbor contacts between monomers of type i and 
j. In this work, we choose €aa = — le, £bb = — 2e, and 
€ab = — 2e, where e is the unit of energy. The large dif- 
ference in the interaction energies e AA and e BB makes 
differences between A and B chains readily observable 
in simulations. The value of e AB corresponds to very 
attractive interactions between unlike monomers and in- 
sures the miscibility of the blends. From the temperature 
T, the unit of energy, e, and Boltzmann's constant, fee, 
a dimensionless temperature T* and its inverse f3 are de- 
fined as T* = k B T/e and ft = l/T*. In this work, we 
present blend simulation results for the fixed tempera- 
ture T* = 10 corresponding to (3 = 0.1. 

For a cubic lattice with L lattice sites on the side, the 
monomer density p is defined as 

N pA N A + N pB N B 
P= jj, , (9) 

where N p a and N pB are the number of chains, while Na 
and N B are the chain lengths for chains of type A and 
B, respectively. The mass fraction ca of component A in 
a blend is given by 

c A = . (10) 

N pA N A + N pB N B v ' 

In this work, we consider only monodisperse melts and 
blends, i.e. N A = N B = N, so that the mass fraction is 
equal to the mole fraction. 

Shaffer established that a monomer density of p = 
0.5a -3 corresponds to a dense melt for the athermal 
system; 5 - We performed Monte Carlo simulations for the 
three monomer densities p = 0.5a -3 , 0.6a -3 , and 0.7a -3 



at (3 = 0.1. For melts and 50/50 blends, we considered 
nine chain lengths between N = 5 and N = 80. For inter- 
mediate concentrations, we focused on chains of length 
N = 20. Simulations were performed on a lattice of size 
L = 20 with periodic boundary conditions applied along 
the three coordinate directions. In order to test for finite 
size effects, we performed simulations for chains of length 
N = 80 in a 50/50 blend on a lattice of size L = 40 and 
found no significant differences in the results. 

Initial configurations were created by randomly placing 
dimers on the lattice and repeatedly reassembling shorter 
chains into longer chains, where the no bond-crossing 
condition was enforced from the start. These initial con- 
figurations were equilibrated at the inverse temperature 
(3 = 0.1 before a trajectory of the simulation consist- 
ing of 10,000 configurations separated by a fixed number 
of Monte Carlo steps, t r , was recorded. The number of 
Monte Carlo steps during production varied between 10 6 
for the shortest chains and 10 8 for the longest chains. 
Typical simulation times for chains of length N — 20 are 
10 7 Monte Carlo steps. In order to improve the statistics 
for blends, where the mass fraction of one of the compo- 
nents is very small (5%), simulation times were extended 
to at least 5 x 10 7 . In all simulations, the equilibra- 
tion time was at least 10% of the production time and 
was sufficiently long for the chains to travel a distance 
corresponding to multiple times the radius of gyration. 
During the simulations, the acceptance rates ^4 ra te for el- 
ementary moves are monitored separately for each of the 
components. 

III. EVALUATION OF SIMULATION RESULTS 

The Monte Carlo simulations are evaluated to yield 
static and dynamic quantities. The average radius of 
gyration, R g , A , of chains of type A, for example, is cal- 
culated from 

N pA 

^-^E^E((^-^) 2 >, (ii) 

pA k=l i<j 

where r^- is the position of monomer i on chain k of 
component A, and where the angular brackets indicate 
the average over configurations in the trajectory. The 
calculation for the radius of gyration of chains of type 
B proceeds in the same way. In order to simplify nota- 
tion, we drop the subscripts indicating the type of chain 
whenever the calculation is identical for both components 
in the blend and there is no danger of confusion. Error 
estimates for static quantities are obtained from block 
averaging^ where we divide the trajectories into ten 
equal blocks. 

In melts and miscible blends, polymer chains of suffi- 
cient length are expected to obey Gaussian statistics. For 
such chains, the average squared end-to-end distance B? e 
is proportional to the number of bonds, B% = b 2 (N — 1), 
where b is the so-called effective bond length. 1 Since the 



4 



squared averages of the radius of gyration and end-to-end 
distance for dense systems are related through R\ = 6R 2 ,, 
we calculate the effective bond length from 



10- 3 



6R 2 J(N-l). 



(12) 



For Gaussian chains, the packing length defined in Eq. © 

may be expressed as p = (pb 2 ) , where p is the 
monomer density defined in Eq. |(9|). This allows us to 
calculate an estimate for the characteristic chain length 
N c from Eq. gj) 



N c = C c (b 3 p) 



(13) 



where the constant of proportionality, C' c , is determined 
from a comparison with the crossover equation (J7J of 
Hessi^a 

In this work, we present results for two mean-squared 
displacement functions. The mean squared displacement 
of the center of mass is obtained from 



1 



.,*(*) 



,,*(o)) 5 



(14) 



where r cm ^ is the position of the center of mass of chain 
k and where the angular brackets indicate the average 
over configurations that are t Monte Carlo steps apart. 
Similarly, the mean squared displacement of the central 
monomer is obtained from 



3i W 



N p 



N/2 



k(t) 



r N/2,k 



(o)) 2 



(15) 



fc=i 



where v N ^ 2 .k represents the position of the central 
monomer of the chain (the results are averaged over both 
innermost monomers for chains with an even number of 
monomers) . 

The self-diffusion coefficients of the components are de- 
termined from the long-time limit of mean-squared dis- 
placement gd of the center of mass 



D 



lim 

t — >oc 



9d(t) 
6t 



(16) 



In Fig. |21 we present simulation results for gd/Qt for both 
components of an A-rich blend (ca = 0.75) of chains of 
length N = 20 at a monomer density of p = 0.7 and an in- 
verse temperature of j3 = 0.1. The results in Fig. |2Jwere 
calculated by two different methods. During the simula- 
tions, a block algorithm adapted from Frenkel and Smit 55 
was used to determine mean squared displacements. The 
algorithm yields values of the displacement functions at 
intervals that increase with increasing time and, thus, 
gives access to a large time range. These values are in- 
dicated by open symbols in Fig. |2 We also calculated 
displacement functions from the trajectories at the time 
intervals t r for times up to half the total run time. These 
results are represented by solid lines in Fig. [3 The graphs 
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FIG. 2: The function ga/Qt, where gd is the mean squared 
displacement and t is the time, as a function of time for an A- 
rich blend with mass fraction ca = 0.75, monomer density p = 
0.7a -3 , inverse temperature = 0.1, and chain length TV = 
20. The symbols indicate results for the two components, A 
(diamonds) and B (circles), obtained with the block algorithm 
during the simulation. The black lines represent results from 
the evaluation of the trajectory, which is limited to times 
larger than the recording interval t r = 10 3 MCs and agrees 
well with the results from the block algorithm. The gray areas 
show the ranges of gd/6t data that were used to determine the 
self-diffusion coefficients of the two components as explained 
in the text. 



show that the values for gd/6t obtained by the two meth- 
ods agree well with each other. The block algorithm is 
very efficient and was used to determine most of the self- 
diffusion coefficients presented in this work. The results 
for gd/Qt decrease with time at short times before they 
level off to a constant value and, finally, become irregular 
at very long times. The decrease of gd/Gt at short times 
indicates the subdiffusive behavior expected for shorter 
chainsi^ For very long times, the results for the dis- 
placement functions are irregular since they represent av- 
erages over few configurations so that individual events, 
like the release of a chain from entanglements, can alter 
the shape of the functions^ The gd/Qt values presented 
in Fig. |3 for the B component are more noisy than those 
for the A component since the blend is rich in A and con- 
tains fewer chains of type B. Self-diffusion coefficients are 
determined from the data in the time range where gd/6t 
is nearly constant; the range starts when gd/Qt closely ap- 
proaches its asymptotic plateau and excludes the largest 
times, where the results for gd/^t become irregular. In 
Fig. [3 we show as gray-shaded areas the data ranges that 
were used to determine self-diffusion coefficients for the 
blend components. For each component, we calculate the 
average of the values of gd/6t in this range, and we also 
fit a function of the form f(t) = Z? + const/6i to the data. 
When the times in the fit range are sufficiently long, the 
second term in f(t) is very small and the values of the 
diffusion coefficients obtained by the two methods agree 
with each other within their statistical uncertainty. In 
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FIG. 3: Simulation results for the squared effective bond 
lengths, b 2 = 6Rg/(N — 1), of polymer chains as a function 
of the number of bonds, N — 1, for a 50/50 blend and for 
the A and B melts at the temperature and monomer density 
indicated in the figure. The upper and lower panels show the 
results for chains of type A and B, respectively. A comparison 
of the chain dimensions in the blends (open symbols) with 
those in the melts (filled symbols) shows that both types of 
chains are expanded in the blend. 



this work, we took care to have simulation runs of suf- 
ficient length to allow for a consistent determination of 
the self-diffusion coefficients. 

In order to obtain estimates for the local friction coef- 
ficients from our simulations, we combine measurements 
of the effective bond length with results for the accep- 
tance rates of elementary moves. According to the Rouse 
model, the segmental mobility [i is related to the shortest 
Rouse relaxation time through 1 

m=A— > ( i7 ) 

where b is the effective bond length and tjv is the relax- 
ation time of a single segment. Since tjv is expected to 
be inversely proportional to the acceptance rate, which 
measures the probability that a monomer is able to com- 
plete an attempted move to a nearest-neighbor site, we 
set in this work 



H = C^b 2 A iatc , (18) 

where A ratc is the acceptance rate. C M is a constant of 
proportionality that is determined from the requirement 
that D/D c = l for N/N c = 1, see Eqs. © and ©. The 
local friction coefficients, finally, are determined from 



IV. RESULTS 
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The effective bond length, which is related to the ra- 
dius of gyration through Eq. I|12|) depends on the struc- 
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FIG. 4: Self diffusion coefficients as a function of chain 
length for all melts and blends considered in this work. The 
filled and open symbols represent simulation results for melts 
and blends, respectively, while diamonds and circles stand for 
chains of type A and B, respectively. The results represent 
three different monomer densities, p = 0.5a -3 , p — 0.6a -3 , 
and p = 0.7a -3 , and show the expected decrease of the self- 
diffusion coefficients with increasing density. In this graph, 
results for the B melts are hidden by the symbols for blend 
results. 



ture as well as the environment of the chains. In Fig.E|we 
present simulation results for the squared effective bond 
length as a function of chain length for A and B chains in 
the melt and in a blend with mass fraction ca = 0.5. We 
note, first of all, that for both melts and blends the bond 
length approaches a constant value for long chain lengths, 
as expected for Gaussian chains, see Eq. I|12l) . Further- 
more, the effective bond length of both A and B chains 
is larger in the blend than in the melt. At all densities, 
the chain expansion increases with increasing dilution, 
i.e. chains are most expanded when they are surrounded 
by chains of the other component. This is due to the 
relatively large attractive interaction energy associated 
with contacts between unlike monomers; the mixed in- 
teraction with energy cab = — 2e is twice as attractive as 
the interaction between A segments, €aa = — le, so that 
A monomers seek out B monomers as neighbors, which 
stretches the chains. 

Self diffusion coefficients give information about the 
large scale dynamics of melts and blends. They are deter- 
mined from Monte Carlo simulations as described in Sec- 
tion lllll In Fig. 0] we present results for the self-diffusion 
coefficients as a function of chain length for all melts and 
blends considered in this work. As expected, the self- 
diffusion coefficients decrease with increasing monomer 
density p and chain length N. The chain-length de- 
pendence of the self-diffusion coefficients in Fig. 0] is 
intermediate between the scaling laws D ~ TV -1 and 
D ~ TV -2 expected from Rouse and reptation theory, re- 
spectively. This is expected since chain-end effects mask 
the Rouse behavi or of short chains in dense systems (see 
e.g. Refs. I58l59|) while fully entangled behavior is not 
usually seen for chains of the relatively short lengths sim- 
ulated herei 42 i 43 i 44 i 45 i 46 i 47 The graph also shows that the 
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FIG. 5: Scaling representation of the self-diffusion coefficients 
for all melts and blends considered in this work. The symbols, 
which are the same as in Fig. 2] represent values of the self- 
diffusion coefficient D divided by D c = p/N c calculated from 
Eqs. jjjj, jIHt i and W'M as a function of the chain length N 
divided by the characteristic chain length N c calculated from 
Eq. il l .'-it . The dashed line represents the crossover function 
of Hessf— see Eq. J7J, while the gray solid lines represent 
the limiting power laws of the Rouse and reptation models, 1 
see Eq. JBJ. 



self-diffusion coefficients for A melts are typically larger 
than those for the blends. (The symbols for the B melts 
are hidden by the symbols for the blends.) The dif- 
ferences between the dynamics of the two components 
and the composition dependence of the self-diffusion co- 
efficients will be discussed in detail for chains of length 
N = 20 below. 

In Fig. we present the self-diffusion coefficient data 
of Fig. 2| in scaling form, D/D c as a function of N/N c , 
where D c is calculated from Eqs. (J5J, l|18|l . and i|13|) 
and N c is calculated from Eq. ifTSjl . In the graph, both 
melt and blend data for different chain lengths and den- 
sities collapse onto a common line. This suggests that 
scaling of dynamic properties with the packing length 
may be applicable to blends as well as melts. It also 
gives us some confidence in our construction of mobil- 
ities from the acceptance rates and the effective bond 
lengths (/i ~ o 2 A ratc ). A comparison with the gray solid 
lines indicating the power laws expected from Rouse and 
reptation theory 1 (Eq. 0) confirms that most of our data 
are in the crossover region between unentangled and fully 
entangled behavior. The dashed line in Fig. [5] represents 
the crossover equation Q of Hess^S*^ which was used 
to determine values for the two constants C M and C c ap- 
pearing in Eqs. I|18|l and Q13[l. There is some uncertainty 
in the results for the constants since more than one com- 
bination of Cfj, and C c values leads to a respectable rep- 
resentation of the data in the crossover region. However, 
a difference in a constant prefactor for the mobilities \x is 
of little concern since we are interested in the variation 
with density and composition rather than the absolute 
values of the mobilities and the corresponding friction 
coefficients £ = l/pf3. 
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FIG. 6: Average mean squared displacement <?i of the in- 
nermost monomer as a function of time for chains of type A 
with three different reduced chain lengths N/N c . The sym- 
bols connected by thin solid lines represent, from top to bot- 
tom, results for chains of length N = 5 in a melt of density 
p = 0.5a~ 3 (N/N c = 0.19), chains of length N = 20 in a melt 
of density p = 0.6a -3 (N/N c = 0.96), and chains of length 
N = 80 in a 50/50 blend of density p = 0.7a -3 (N/N c = 6.2). 
The dashed lines indicate the power laws expected for gi(t) 
from reptation theory^ 



To help interpret our values for the characteristic chain 
lengths N c we present in Fig. [5] average mean-squared 
displacements of the central monomer, gi(t), for three 
systems corresponding to the smallest, the largest and 
an intermediate value of N/N c in Fig. [S] The sym- 
bols connected by solid lines represent values calculated 
from Eq. I|15fl with the block algorithm. They repre- 
sent results for an A-melt of chains of length N = 5 
at density p = 0.5a" 3 with N/N c = 0.19 (top), an A- 
melt of chains of length N — 20 at density p — 0.6a -3 
with N/N c = 0.96 (center), and the A component in 
a 50/50 blend of chains of length N = 80 at density 
p = 0.7a~ 3 with N/N c = 6.2 (bottom). All three exam- 
ples show the expected diffusive behavior (<?i ~ t) at long 
times. The behavior at short times is subdiffusive but not 
quite Rouse- like (gi ~ t 1 / 2 ), as is generally observed in 
simulations^*^ The largest gi(t) values correspond to 
the smallest reduced chain length, N/N c — 0.19, and il- 
lustrate unentangled behavior; the slope of g\(t) in this 
double logarithmic plot increases monotonically from the 
subdiffusive region at short times until it reaches a value 
of unity in the diffusive region, which starts around 10 3 
MCs for this system. The smallest gi(t) values corre- 
spond to the largest reduced chain length, N/N c = 6.2. 
In this second subdiffusive region at intermediate 

times is clearly visible. The slope in this region has a 
value of about 0.37 which is not as small as the value 
of 1/4 expected for fully entangled chains from repta- 
tion theoryi This implies that even the longest chains in 
our simulations are not fully entangled. The line of in- 
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FIG. 7: Self-diffusion coefficients of melts and blends as a 
function of blend composition for chains of length N — 20 
and monomer densities p = 0.5j4 -3 (top), p = 0.6A -3 (cen- 
ter), and p — 0.7j4 -3 (bottom). Results for A melts and the 
A component in blends are represented by open diamonds, 
while those for B melts and the B component in blends are 
represented by filled circles; ca denotes the mass fraction of 
component A. 



FIG. 8: Scaling representation of the self-diffusion coeffi- 
cients for melts and blends with chain length = 20 at the 
monomer densities p — 0.5a -3 , p — 0.6a -3 , and p = 0.7a -3 . 
The open diamonds and filled circles correspond to results for 
the A and B component, respectively. For each density, the 
results for the pure melts correspond to the lowest reduced 
chain lengths N/N c and the largest values of the scaled self- 
diffusion coefficients D/D c . Conversely, the largest values of 
N/N c and the smallest values of D/D c for each density cor- 
respond to the minority component in a dilute blend. 



termediate g\ (t) values in Fig. H3 finally, corresponds to 
a melt with chains near the characteristic chain length 
(N/N c = 0.96). In this case, the slope of gi(t) de- 
creases only slightly before it rises to the diffusive value 
of unity. This suggests that the chains are just starting to 
feel entanglement effects for reduced chain lengths near 
N/N c ~ 1. 

In order to investigate the composition dependence of 
the dynamics, we focus in the following on one chain 
length, TV = 20, and investigate a broad range of com- 
positions ranging from the pure melts (ca = 0, ca = 1) 
to the dilute limits (ca — 0.05, ca — 0.95) with three 
intermediate compositions (ca — 0.25, ca — 0.5, and 
ca = 0.75), for the three monomer densities, p = 0.5a -3 , 
p = 0.6a -3 , and p = 0.7a -3 and reduced inverse tem- 
perature /3 = 0.1. In Fig. we present the simulation 
results for the self-diffusion coefficients of the N = 20 
melts and blends as a function of composition for the 
three monomer densities p — 0.5a -3 , p = 0.6a -3 , and 
p = 0.7a -3 . These results show that, in general, the 
composition dependence of the self-diffusion coefficients 
decreases with increasing density and is stronger for the 
A component than for the B component. The trends for 
the self-diffusion coefficients are the same at each density; 
as the mass fraction ca of the A component increases, the 
self-diffusion coefficients of A chains increase while those 
of B chains decrease. Hence, for both types of chains, 
blending tends to reduce the values of the self-diffusion 
coefficients from the melt values (the self-diffusion coeffi- 
cients of the B components are largest in B melts and de- 
crease as A component is added, while the D values of A 
chains are largest in the A melt and decrease as B compo- 
nent is added to the blends) . For most blends considered 



here, the self-diffusion coefficients of the A component are 
larger than those of the B component. This is true for all 
blends of the lowest density considered here, p = 0.5a -3 . 
However, for the density p — 0.6a -3 , the values of the 
self-diffusion coefficients of both components are compa- 
rable for blends with small A content. For the highest 
density considered in our simulations, p = 0.7a -3 , finally, 
the self-diffusion coefficient of the A component is even 
smaller than that of the B component for the smallest 
blend value of ca ■ 

In Fig. [SJ we present simulation results of the self- 
diffusion coefficients in scaled form, D/D c as a function 
of N/N c . The data are a subset of those presented in 
Fig. [5] and belong to the crossover region near N/N c = 1. 
In the double-logarithmic graph of Fig. HO the reduced 
self-diffusion coefficient data follow closely a straight line 
with a slope of about -1.7. The reduced diffusion coeffi- 
cients of the minority components in dilute blends appear 
to lie somewhat above this line; however, this may not 
be significant considering the uncertainty in the data. 
Hence, the data in Fig. |S] may be approximated by 

D/D c ~ const, x (N/N c y v with v ~ 1.7, (20) 

which implies 

D ~ const, x p,N"~ l N~ v , (21) 

where D c = p,/N c of Eq. (JSJ has been used. The results 
presented in Fig.|S|suggest that the scaling with a charac- 
teristic chain length N c derived from the packing length 
and a characteristic diffusion coefficient D c — pi/N c con- 
structed with the mobility p ~ b 2 A [ate captures both 
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FIG. 9: Scaling chain length N c as a function of blend com- 
position for melts and blends with chain length N = 20 at 
the indicated densities. The symbols connected by solid lines 
represent values calculated from Eq. I13H with Eq. 112> for 
components A (open diamonds) and B (filled circles). 

density and composition dependence of the self diffusion 
coefficients for these blends. 

In order to gain insight into the origin of the vari- 
ation of the self-diffusion coefficients with density and 
blend composition, we present in Figs. El El and ^2 re- 
sults for the characteristic chain lengths, the acceptance 
rates, and the mobilities of the melts and blends with 
N = 20. The results for the scaling chain lengths N c 
m Fig. show that N c increases with increasing density. 
This suggests a decrease of the entanglement length with 
increasing density, which is expected since the distance 
between monomers from different chains decreases with 
increasing density 4i The N c values in Fig. ED also show 
that blending decreases the characteristic chain length 
for the blends considered here; for both types of chain, 
the values of N c decrease with increasing concentration 
of the other component at a given density. This com- 
position dependence of N c is a result of the variation of 
the effective bond length b with composition. The results 
for b 2 in Fig. 1 1 21 show that both A and B chains are ex- 
panded in a blend compared to the melt. Since, according 
to Eq. (fH5|) , a larger effective bond length corresponds to 
a smaller characteristic chain length, the values of N c 
decrease upon blending from their melt values. 

The average acceptance rates for elementary moves of 
our Monte Carlo simulations presented in Fig. ^3 mea- 
sure the average probability for a monomer to complete 
an attempted move to a nearest-neighbor site. In the lat- 
tice model for polymer blends employed in this work, this 
probability depends, first of all, on the local bond struc- 
ture since only three different bond lengths are allowed. 
It also depends strongly on the local density since the 
monomer cannot move to a site that is already occupied 
by another monomer. The composition of the blends af- 
fects the acceptance rate through the change in internal 
energy associated with a move. The internal energy of 
the system depends on the numbers (Nij) and the inter- 



18 I — L 




6 — i 1 1 1 1 1 — 

0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 10: Average acceptance rates for elementary moves as a 
function of blend composition for melts and blends with chain 
length N = 20 at the indicated densities. The open diamonds 
and filled circles connected by solid lines represent simulation 
results for components A and B, respectively. 



action energies (ey, i,j <G {A, B}) of contacts between 
non-bonded nearest neighbors, see Eq. |JS}. According to 
the Metropolis criterion, a move is always accepted when 
the new configuration has a lower internal energy than 
the original one. If the new internal energy is higher, the 
move is accepted only with a probability corresponding 
to the Boltzmann factor of the energy difference between 
old and new configurations. This implies that the ac- 
ceptance rate for moves is small, when the new site has 
a smaller number of occupied nearest neighbor sites or 
when the interactions with the new neighbors are less 
attractive than those with the old. Finally, a small frac- 
tion of moves is prohibited because it would lead to the 
crossing of bonds. The results for the acceptance rates 
in Fig. lfTU|) show the expected strong dependence on the 
monomer density. In order to understand the composi- 
tion dependence, we start by considering blends, where 
the A component is dilute, i.e. near ca = 0. In this case, 
the A chains are surrounded by B chains. Since the in- 
teraction energies for AB and BB contacts are identical, 
(-ab = £bb = — 2e, we expect similar acceptance rates for 
the A and B component at a given density. This is in- 
deed what the data in Fig. 1 101 show. For all compositions 
and densities, the acceptance rates for the A component 
are higher than those of the B component. Furthermore, 
for a given density, as A component is added to a blend, 
the mobility of A monomers increases while that of B 
monomers decreases. There are two factors, both related 
to the energetics, that contribute to this. First of all, 
since A-A contacts are less attractive than A-B contacts 
{caa = — le while €ab = — 2e), A monomers become in- 
creasingly more mobile as A component is added to a 
blend. Secondly, the local density near A monomers is 
somewhat lower than that near B monomers leading to a 
larger mobility of A monomers. The reason for the local 
density variation is the difference in the average energy 
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FIG. 11: Mobility pas a function of blend composition for 
melts and blends with chain length N = 20 at the indicated 
densities. The symbols connected by solid lines represent val- 
ues calculated from Eq. 1181 with Eq. (I I 12 P for components A 
(open diamonds) and B (filled circles). 



penalty for making a contact with a void instead of an- 
other monomer. For B monomers, this penalty is always 
2e. For A monomers, on the other hand, it varies with 
composition between 2e near ca — and le near ca = 1- 
The effect becomes smaller with increasing monomer den- 
sity since the number of voids decreases. 

In Fig. ^] we present results for the composition de- 
pendence of the mobilities, calculated from Eq. (|18fl with 
Eq. i|12|) . for the melts and blends with chain length 
N = 20. The mobilities are proportional to both the 
acceptance rates and the squared effective bond lengths, 
p r~j b 2 A Yat o- It is not surprising that they show some 
of the characteristics of the acceptance rates discussed 
above. Just like the acceptance rate, the mobility de- 
creases with increasing density. Furthermore, chains of 
type A always have a larger mobility than chains of type 
B, which makes A the faster component in the blend at 
all compositions. However, the composition dependence 
of the effective bond length partially compensates that 
of the acceptance rate and makes the mobilities less de- 
pendent on composition than the acceptance rates. As A 
component is added to a blend, the mobility of A chains 
increases, while that of B chains decreases only slightly. 
In the limit where the A component is dilute, the mo- 
bilities of the two components remain distinct, while the 
acceptance rates of the two components converge. Con- 
versely, in the limit of dilute B in A, the relative differ- 
ences between mobilities of the components are smaller 
than those of the acceptance rates. 

The local friction coefficients £ of the blend compo- 
nents are inversely proportional to the mobilities, £ — 
l/f3fi (Eq. CHJ)- I n the left panel of Fig. El w e present 
values for the local friction coefficients calculated from 
the mobility results for melts and blends of density 
p = 0.6a -3 presented in Fig. ^] Since A is the faster 
component, its friction coefficients are smaller than those 
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FIG. 12: Local friction coefficients as a function of composi- 
tion for melts and blends with chain length N = 20. The sym- 
bols in the left panel represent values £ = correspond- 
ing to the mobilities p at the constant density p = 0.6a~ 3 
in Fig. 1111 The symbols in the right panel represent friction 
coefficient values interpolated from the data in Fig. 1131 to a 



constant pressure of P ■ 
the symbols. 



2e/a . The solid lines only connect 



of the B component. As the mass fraction of A chains 
increases, the friction coefficients of the A component de- 
crease (slightly) while those of the B component increase. 
Hence, when A component is added to a blend at con- 
stant density, A segments speed up (somewhat) while B 
segments slow down. This is the case for all densities con- 
sidered here and contrasts with the experimental obser- 
vation that both components speed up when the fraction 
of faster component in a blend is increased. The rea- 
son for this apparent discrepancy is that we have varied 
the composition at constant density p, whereas in ex- 
periments the composition is typically varied at constant 
pressure. In order to arrive at friction coefficient values at 
constant pressure, we first consider their pressure depen- 
dence and then interpolate our data to a given pressure. 
In this work, we estimate the pressure using a random 
mixing approximation for the internal energy and the 
Flory-Huggins expression for the entropic contribution. 60 
This is expected to yield at least qualitatively correct re- 
sults since our simulations are carried out at high density, 
where the Flory-Huggins theory has been found to give 
a reasonable representation of the pressure of athermal 
lattice chainSf£ii§£ and since we work at a constant, high 
temperature where the random mixing approximation is 
expected to be adequate. In Fig.^Jwe present results for 
the friction coefficients of melts and blends with N = 20 
as a function of pressure. As expected, the friction coeffi- 
cients increase with increasing pressure for each composi- 
tion. In this representation it becomes apparent that the 
friction coefficients of both components decrease when 
the mass fraction of A is increased at constant pressure. 
To illustrate this further, we have estimated values for 
friction coefficients at the pressure P — 2e/a 3 by inter- 
polation of the results in Fig. These values are pre- 
sented as a function of mass fraction in the right panel of 
Fig. 1121 In qualitative agreement with experiment, the 
results for the friction coefficients show that the dynam- 
ics of both components speed up when the fraction of the 
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FIG. 13: Local friction coefficients as a function of pressure 
for melts and blends with chain length N = 20. The sym- 
bols represent values of £ = 1//3/J, corresponding to the mo- 
bilities p in Fig. 1111 The pressure has been estimated from 
Flory-Huggins theory p— The dashed and dashed-dotted lines 
connect values for the B and A components, respectively, in 
blends of a given composition; melt values are connected by 
solid lines. The results shown here correspond to mass frac- 
tions c A = 0, 0.05, 0.25, 0.5, 0.75, 0.95, and 1.0. The friction 
coefficients for the B component in the melt and in the blend 
with ca = 0.05 are too close to be distinguished in this rep- 
resentation. 

fast component (A) is increased at constant pressure. 

V. DISCUSSION 

In this work, we have presented Monte Carlo simula- 
tion results for polymer melts and blends with a range 
of densities, compositions and chain lengths. In order to 
represent the two components of a binary blend, we have 
generalized Shaffer's bond fluctuation model for ather- 
mal melts£i*£2iS by introducing attractive interactions 
between occupied non-bonded nearest neighbor sites on 
the lattice. The interactions parameters employed in this 
work are caa — — le and €bb = — 2e for like interactions 
between monomers of type A and B, respectively, and 
eab = — 2e for interactions between unlike monomers. 
All simulations were carried out a fixed temperature of 
T = lOe/fce corresponding to an inverse reduced tem- 
perature of (3 = 0.1. Our results for the effective bond 
length presented in Fig.[5]show that the chain dimensions 
of both types of chains increase upon blending. This ef- 
fect increases with increasing dilution and is due to the 
large attractive interaction energy associated with con- 
tacts between unlike monomers. 

Self-diffusion coefficients were determined from the 
mean square displacement of the center of mass of the 
chains. The simulation results presented in Fig. 0] 
show that the chain-length dependence of the self- 
diffusion coefficients falls between the power laws ex- 
pected from Rouse and reptation theory^ which is ex- 
pected for the range of chain lengths considered in this 



work. 42,43, ■ 5-46 ' In order to investigate the scaling 
behavior of the self-diffusion coefficients, characteris- 
tic chain lengths N c were calculated from the packing 
lengths of the chains. Segmental mobilities p were deter- 
mined from the acceptance rates and the effective bond 
lengths and combined with the values for N c to yield 
characteristic self-diffusion coefficients D c — p/N c . The 
results for the reduced diffusion coefficients D/D c as a 
function of the reduced chain length N/N c presented in 
Figs. [S] and Figs. [S] show the data for melts and blends 
with different chain lengths, densities, and compositions 
to collapse onto a common line. This suggests that the 
scaling of dynamic properties with the packing length, 
which has been observed for polymer melts^SiiSL^&SLiSdi 
may be applicable to blends as well. 

The composition dependence of the dynamic proper- 
ties was studied in detail for chains of length N = 20 in 
melts and blends of three different densities. The results 
for the self diffusion coefficients D presented in Fig. 
show the expected variation with density. At each den- 
sity and for both types of chains, blending reduces the 
self-diffusion coefficient from its melt value. This may 
be attributed to the composition dependence of the mo- 
bility p and the characteristic chain length N c , as fol- 
lows. The scaling representation of Fig. [S] for chains of 
length N = 20 shows that the reduced self-diffusion co- 
efficient follow approximately a power law with exponent 
v = 1.7. This implies that D depends on the mobility 
p and the characteristic chain length N c approximately 
as D = const, x pN"~ 1 N~ u , see Eq. ffiQ. Since the 
values of both the mobility and the characteristic chain 
length decrease upon blending, so do the values of the 
self-diffusion coefficients. 

The variation of the self-diffusion coefficients with com- 
position has another interesting aspect. For the high- 
est density, p = 0.7a -3 the difference between the self- 
diffusion coefficients of the two components changes sign 
from Da — Db < to Da — Db > at a mass fraction 
of A of about ca — 0.2. This is surprising since it im- 
plies that, for p = 0.7a -3 , A chains move more slowly 
than B chains at low A concentrations but faster than B 
chains at higher A concentrations. Since changes in both 
mobility and characteristic chain length affect the value 
of the diffusion coefficient we consider them in turn. Re- 
sults for the segmental mobility presented in Fig.llllshow 
that the local mobility of the A chains in a blend is always 
larger than that of the B chains. This implies that, on 
the length scale of an effective bond, the A component is 
always the fast component in the blends considered here. 
Hence, the segmental mobility cannot be responsible for 
the sign change of Da — Db- Results for the characteris- 
tic chain lengths N c in Fig.[§|show that, for all densities, 
the differences between the characteristic chain lengths 
change sign from N c a — N c b < to N c a — N c b > at 
intermediate mass fractions. This, together with the rel- 
atively weak composition dependence of the mobility at 
p = 0.7a" 3 , gives rise to the sign change of Da — Db- In 
physical terms; for the highest density considered here 
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and at the lowest concentration of A chains, the char- 
acteristic chain length for A chains is so much smaller 
than that of B chains, that entanglement effects lead to 
A chains diffusing more slowly than B chains even though 
the segmental mobility of the A chains is larger than that 
of the B chains. 

Values for the local friction coefficients ( have been 
determined from the results for the segmental mobilities 
fi. The pressure dependence of the friction coefficients 
is shown in Fig. 1131 where the values for the pressure 
are estimated from the Flory-Huggins theory 1S . The re- 
sults show the expected increase in the friction coeffi- 
cients with increasing pressure. The pressure dependence 
is slightly larger for the slower (B) component than for 
the faster (A) component, which agrees qualitatively with 
results from our calculations for polyolefm blends^* In 
Fig. El we contrast two sets of results for the composition 
dependence of the local friction coefficients. The results 
in the left panel, where the composition of the blends is 
varied at constant density, as is typical for simulations, 
are qualitatively different from those in the right panel, 
where the pressure is kept approximately constant, as is 
typical for experiments. The results for the constant- 
pressure composition dependence of the friction coeffi- 
cients agree qualitatively with experimental observations 
on miscible blends; the friction coefficients of both com- 
ponents decrease with increasing mass fraction of the fast 
component (A). 

In this work, we calculate the local friction coefficients 
from the mobility which describes the dynamics on 
the scale of the effective bond length. For polymer melts 



near the glass transition^ and for miscible polyolefin 
blends^ the Kuhn length has been identified as the 
size of the moving segment relevant to local dynamics. 
For the blend model considered here, the effective bond 
length and the Kuhn length are comparable in size, which 
makes it difficult to distinguish between these choices for 
the size of the moving segments. However, we hope that 
a current investigation of the dynamics of athermal melts 
for a larger set of densities will give us some insight into 
this question^ 

It is encouraging that these first simulation results for 
our lattice model for polymer blends reproduce trends ob- 
served in experimental work on polymer blends. The dif- 
ferences between the friction coefficients of the two melts 
and between the melts and the blends, however, are fairly 
small. The reason is that chains of type A and B differ 
only in the values of the interaction parameters and that 
the simulations are carried out at high temperature. In 
real polymer blends, the blend components typically dif- 
fer in chain stiffness, which is known to have a large effect 
on the dynamics iSSi 3 ^^! 3 . In order to simulate more re- 
alistic systems and to increase the difference between the 
chains of the blend components, we are currently extend- 
ing our model to include a differences in chain stiffness. 
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